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I. Summary 

When a simple model for the relationship between the density- temperature 
fluctuation correlation and mean values is used, we determine that the rate of 
change of turbulent intensity can influence directly the accretion rate of droplets. 
We find, experimentally, that the droplet particle size is both temperature and 
Reynolds number dependent. We also find that the rate of droplet growth has an 
important dependence on Reynolds number, in a manner stronger than its sim- 
ple temperature dependence, thereby suggesting a previously unsuspected role 
for turbulence in the control of condensate accretion. 

II. Background 

Heterogeneous nucleation and the subsequent droplet and crystal growth 
play a critical role in the contribution from PCS and vapor trails to atmospheric 
environmental concerns and in the development of ice crystals on environmen- 
tally exposed aerofoils. Yet, the best growth rate theories 1 ' 2 ' 3 ' 4 are essentially ad 
hoc and/or empirical with very weak experimental confirmation 5 ' 6 of underlying 
physical principles. In a recent work, 7 U. De Silva , A. Gardner and J. A. Johnson 
HI have shown experimentally that the droplet size increases at fixed Reynolds 
number when the temperature decreases. They have also shown, independent 
of temperature, a pronounced dependence of droplet size on the Reynolds num- 


1 Lothe ,J.,and Pound,G.M., "Reconsideration of Nucleation Theory", Journal of Chemical Physics, 
Vol.36, No.8, 1962, pp.2080-2085. 

2 Lothe,J.,and Pound,G.M., "Statistical Mechanics of Nucleation", in Nucleation edited by A. C. 
Zettlemoyer, Marcel Dekker, New York, 1969, Chap. 3. 

3 Dunning, W. J., "General and Theoretical Introduction," in Nucleation edited by A. C. 
Zettlemoyer, Marcel Dekker, New York, 1969, Chap.l. 

4 Wegner,P.P. "Gas Dynamics of Expansion Flows with Condensation and Homogenous 
Nucleation of Water Vapor" Nonequilibrium Flows, Pt.l, Vol.l, edited by P. P. Wegner, Marcel 
Dekker , New York, 1970, Chap. 4. 

5 Brian, A. B.,Zuev, A. P. , and Khmelevskhi , A. N. /'The effect of Condensation on the 
Parameters of the Flow behind Shock Wave in Water Vapor" Zhumal Tekhnicheskoi Fiziki, 
Vol.62, July 1992, pp.765-770. 

6 Schnerr, G. H. and Bohing, R. "Compressible Turbulent Boundary Layers with Heat Addition by 
Homogenous Condensation", ALAA Journal Vol.30, No.30, 1992, pp. 1284-1289. 

7 U. De Silva, A. Gardner, J. A. Johnson HI, ALAA Journal, 33, 1995, pp. 368-370. 
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ber, specifically suggesting an explicit sensitivity of droplet size to turbulence. In 
that work, the measurements relied on Mie scattering at a single port at ten dif- 
ferent Reynolds numbers and at three different temperatures. Here we will pro- 
vide a theoretical introduction to the problem of droplet size growth rate sensi- 
tivity and then report on our preliminary determinations of a sensitivity in 
droplet size growth rate (by using two ports) to both temperature and turbulence 
explicitly. * 

III. A Theoretical Model for Turbulent Condensate Accretion 
Considerable interest exists in the accretion rate for condensates in 


nonequilibrium flow with icing. 8 Turbulence is thought to play an important 
role in such flows. It has already been experimentally determined, as cited 
above, that turbulence influences the sizes of droplets in the heterogeneous nu- 
cleation of supersaturated vapors. This section addresses the issue of the possi- 
ble influence of turbulence on the accretion rate of droplets. 


According to a nonequilibrium model developed earlier 9 the droplet growth 


rate for flow with changing local pressure is given by: 

= -3BN — log — 
dt n dt 6 P 


with B= (2a)/[p c T(R/p v )] where <r = surface tension coefficient; p c = condensate 
density; and p v = molecular weight of vapor. N n is a correction coefficient such 
that 0<N n <4 and N n =[(n-3) n ' 3 V[(n-4) n * 2 ] for n>4. is obtained through the 
Clausius-Clapeyron equation. For adiabatic flow of a compressible ideal gas the 


mean value of the square of the local speed of sound is defined as 



( 2 ) 


®T. Cebeci, AIAA Journal, 22/ 1995, pp. 1351-2 

9 R. Hazoume and J. Chabi Orou, "Pressure Sensitivity in the Droplet Growth Rate in 


Nonequilibrium Condensation," in Mathematical and Experimental Developments for 
Noneouilibrium Statistical Systems. Institut de Mathematiques et de Sciences Physiques, 


Universite Nationale du Benin, May 25-27, 1994, Contonu, Benin (to be published) 
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where: y = C p /C v is the ratio of the specific heats; g is the mean value of the 

condensate mass fraction of water; P is the mean pressure function; sis the 

mean entropy; and p is the mean density function. 

We can readily write from (2) that: 

p p 

— = cnst and — = cnst (3) 

P (^(TT) 


Equations (2) and (3) were written assuming there are fluctuations of the princi- 
ple physical properties around their mean values as, e.g., g- g + g'with ^ = 0. 
However g ' 2 0. This means that we can consider P — P + P ' , p = p + p' , g 

= g + g' , and T = T + T' . The equation of state 


P = P( 


1 ~(Oq 
Pi 


+ 


<- Qp-g 

p v 


)RT 


( 4 ) 


becomes 


P + P' = (p + p')^—^+ C ^ 

Pi P v 


8-8 


)R(T + V) 


( 5 ) 


with 6) 0 =the initial humidity, p t =the molecular mass of the carrier gas, and 

R=the ideal gas constant. Averaging equation (4) yields: 

p = (i + — ° --£)/?(fp + pT 7 ) (6) 

Pi Pv 

where a new function, p'T' , has appeared; this is the density-temperature corre- 
lation function involving the fluctuations of density and temperature. 

Now we make the important simple assumption 

~p 7 T' = -XpT. (7) 

The equation of state then becomes 

!■ = /(«)(! -A)fff = — (8) 

p r 


where 

l-CQp | CQ 0 -g 

Pi Pv 



( 9 ) 
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We must also make it clear that the physical properties of interest are functions of 
space and time. That is, 

f = f(x i ,t),P = P(x i ,t),p = p(x i ,t),g=g(x i ,t),X = X(x i ,t) (10) 
where i = 1,2,3, xi is the i component of the space variable, and t is the time. To 
use equation (10), the simplest adequate model for A(xi,t) would be 
— = exp— ft. With this model, the density temperature correlation func- 
tion depends only on time and obeys a decay law which is exponential. Finally, 
it is then usual (with U and A as velocity and length parameters respectively) to 
express all these relations in a nondimensional form: V[ =u j/U; y 2 -x,7 A; and 
T=tU/A. 

We can now include the influence of fluctuations in temperature (as well as 

the implicit fluctuations in P and p) in the mean square of the local speed of 

sound a 2 (y 7 ,r ) by explicitly requiring ^ ^ = — . This form is reconciled with 

P(y 2 ^) 7 


equation (8) above by 


b 2 a 2 


* fL_ = ^_(l_A0. 2 ,T)) 

7 7 

y 

We have shown in equation (3) that Pa{b 2 ) (y ~ x) . Thus 


= SJ2J 3T X *Tt U, ^- X ' ) 


= - Z -l-fa 2 (y 1 ,r) f ——) 

y-1 oT exp/f + 1 

in which /is the real positive constant introduced in the definition of X. We fi- 
nally obtain from equation (1) 


6 cr (n- 3) n 3 7 

/ 

da 2 (y 2 ,T) 

ft( _R )f( „-4,- (7-1, 

exp fr + 1 

dr 
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The interpretation of equation (13) is straightforward. (i)At small values of 
x and constant turbulent intensity, ( da 2 (y 2 ,z)/dz=0), the rate of change in droplet 
size is determined entirely by the density correlation exponential factor/. (ii)At 
large values of % and constant turbulent intensity, (da 2 (y 2 ,z)/ dz= 0), the rate of 
change in droplet size has no dependence on turbulence at all. (iii)If the turbu- 
lent intensity is not constant, i.e., (da 2 (y 2 ,z)/ dz^Q), then the rate of change in 
droplet size is sensitive to the turbulence at all times. The case da 2 (y 2 ,z)/ dz*0 
requires that a measure of the actual evolution in the strength of turbulent fluc- 
tuations be determined and that a model for determining the connection between 
the mean square local speed of sound and the evolution of turbulence be deter- 
mined. It has already been shown 2 that, at fixed thermodynamic conditions, the 
size of a droplet in condensing flow is influenced by the value of the Reynolds 
number and, by inference, the strength of local turbulence. Equation (13) above 
sets conditions for the sensitivity of the condensate accretion to turbulence, based 
on basic physics principles, which can now be tested by experiments. 

IV. Experimental Results 

To realize the purpose of droplet rate of growth measurement, we looked 
at the scatter signals from two ports 30 centimeters apart along the shock tubes 
axis each other for ten different Reynolds numbers and five different tempera- 
tures. A pressure driven shock tube was used consisting of a driver chamber of 
60x5.2 in. and a two part driven chamber with a test section of 60x5.2 in. and an 
extension tube of 76x5.2 in. Aluminized mylar sheets of two different thicknesses 
were used as the material for the diaphragm which separates the driven chamber 
from the driver chamber. Compressed air is supplied to the driver chamber from 
a reservoir, where the relative humidity of the air could be controlled by setting 
the temperature in the reservoir at 70°F. 
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A heating tape is wrapped 

around the driver chamber, and a 

thermocouple probe is used in the 

driver chamber to measure the initial 

temperature. The test section has three 

stations for optical diagnostics, two of 

which consist of four perpendicular 

glass ports, the other station has two 

in-line glass ports. An Argon laser is 

used for the right angle scattering through two of the three optical stations. Two 

pressure transducers are mounted on either side of the optical station . The static 

pressure and the right angle scatter signals are recorded on the four channel 500- 

MHz Tektronix digital oscilloscope. All of the signals are stored on disks. The 

optics for receiving the scattered signal from the droplets include a natural 

density filter; the received signals are transmitted to the photomultiplier tubes 

«> 

using the fiber optics and also stored on the oscilloscope just mentioned. 

Appropriate analytical and calibration procedures are used so as to provide 
continuous recordings of droplet size and density at the measuring stations. Off- 
line reconstructions are performed, using appropriate computer software and 
hardware, so as to provide histories in the droplet frame of reference using the 
velocity calculated from the assumptions of one-dimensional shock tube flow 
and confirmed by direct LDV. From these results, we obtain droplet growth rate 
data and a droplet growth rate model which will determine the roles of tempera- 
ture and suspensions in the evolution of droplets. 





0x10° 5x10 7 1x10* 2x10® 2x10 8 3x10 8 
Reynolds numbers 

Figure 2 


We know that, in a pressure driven shock tube, at a shock wave velocity of 
Mach number of 1.8 the contact surface is fully turbulent. 10 The shock wave ve- 
locity is fixed (M=1.8), by maintaining 
the pressure ratio (P4/P1) constant. 

15 

We changed Reynolds numbers by ■ 

1 10 : ______ 

changing Pi and P4 at (P4/P1) _ “ 

£ 5 ; . 80 °F 

constant. After each run, we f /rrAt- . 8 s°f 

I 0 - *W ♦ 1 85 F 

measured the shock wave velocity “■ , s : Sr * 90 ° F 

I : / n ioo°f 

using the quartz pressure transducers. 8 , 10 : J 

The scatter signals from the two .is : 

u 7 1 f i 1 | 1 1 1 1 | 1 1 1 1 | m 1 1 I 1 1 » 1 | 

0x10° 5x10 7 1x10* 2x10* 2x10 e 3x10 8 

photomultipliers were converted to Reynolds numbers 

1 Figure 2 

the particle size by using Raleigh's 

scattering theory. 11 The intensity (in voltage) of the scattered wave at the 

distance r from the origin of axes to the particle is given by Rayleigh formula: 

2 

4 6 2 

2~~4 277 ( 1+cose ) « 

r X 4 m +l 

where a is the particle 's radius ,m is a complex relative refraction index, X is 
the wave length in the medium. For simplicity in making comparisons, we used 

the relative particle size respect to the first particle. Hence the relative particle 
size is given by = 

The influence of the turbulence on the particle size is shown on the Fig.2, 
where we plotted the relative particle size change against Reynolds number. The 
relative particle size change is given by: (— — ). where a is the particle radius, 

a io 

the first index refers to the port number (one or two) and the second one is the 

10 Johnson, J. A. ,Santiago J.and Jones ,W. R. "Driver Gas Flow with Fluctuations", Journal of 
Physics D : Applied Physics,Vol.l3 , 1980 pp.1413-1425. 

■^Kerker,M." The Scattering Of Light." Academic Press, New York,1969, Chap.13 


4 6 

2 

87T a 

m -1 

V 

r X 

tn 2 + 1 
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index for the Reynolds number. The direct explicit connection between turbu- 
lence and Reynolds number and the role of this connection in nonequilibrium 
phenomena have already been shown .12 Fig. 2 shows that the particle size in- 
creases as the temperature decreases at fixed Reynolds number. The data in Fig.2 
are results from the first port. The same results are available from the second 


port. 

We have determined the relative particle size change per unit time as func- 
tion of temperature and Reynolds number. The relative particle size change per 

, a 2j~ a u , 

unit time is given by *- > where t is the time required for the contact surface 

to travel from the port one to the port two. With these measurements, we have 
found a new clear dependence of droplet size on the strength of the turbulence. 
As Fig, 2 above shows, increasing Reynolds number produces an increase in the 
relative sizes of the droplets in a manner which changes with temperature and 
changes with Reynolds number according to a standard strong coherence model. 
Furthermore, when the- rate of droplet growth is measured by observing the 
change in droplet size at two locations, we find that the rate of change of droplet 


• 

T3 a 173K 

o 

T3 a 176K 

■ 

73 a 178K 

▲ 

T3 a 180K 

T 

T3 a 183K 


i t % 


9.0xlb" 7 i.o’xlb ’ 8 1.1x1 6 8 1 ixio Bl.3xf 
Reynolds No. pjg g 


We show a temperature sensitivity in 
these results also. We expect, with this 
motivation and these results, to extend 
our investigations to low temperature 
regimes; we are particularly interested 
in the pressure and temperature condi- 


12 Johnson, J. A. and Ramaiah , R. I. L., "Reduced Molecular Chaos and Flow Instability " 
Stability in the Mechanics of Continua , edited by F. H. Schroeder, Spring-Verlag , Berlin , 1982 
,pp. 318-329. 
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tions appropriate to sublimation. Under these conditions we expect to determine 
the role of both turbulence and temperature in the evolution of ice crystals, ap- 
propriate to high altitude jet trails and to the onset of icing on the wings of air- 
planes. 
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